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A supersonic beam of metastable neon atoms has been decelerated by exploiting the interaction 
between the magnetic moment of the atoms and time-dependent inhomogeneous magnetic fields in 
a multistage Zeeman decelerator. Using 91 deceleration solenoids, the atoms were decelerated from 
an initial velocity of 580 m/s to final velocities as low as 105 m/s, corresponding to a removal of 
more than 95 % of their initial kinetic energy. The phase-space distribution of the cold, decelerated 
atoms was characterized by time-of-fiight and imaging measurements, from which a temperature of 
10 mK was obtained in the moving frame of the decelerated sample. In combination with particle- 
trajectory simulations, these measurements allowed the phase-space acceptance of the decelerator 
to be quantified. The degree of isotope separation that can be achieved by multistage Zeeman 
deceleration was also studied by performing experiments with pulse sequences generated for ^''Ne 
and ^^Ne. 



I. INTRODUCTION 

The development of methods to produce atomic and 
molecular samples with precisely controlled translational 
and internal degrees of freedom has received consider- 
able attention in recent years, with particular emphasis 
on the generation of cold (T = 1 mK — IK) and ultra- 
cold (T < 1 mK) samples. This interest in the study of 
cold molecules is motivated by applications in a broad 
range of disciplines [TJ |2], including the study of colli- 
sion processes and chemical reactions in a temperature 
regime where quantum-mechanical phenomena may dom- 
inate the reaction dynamics, high-precision spectroscopy, 
quantum information processing and quantum simula- 
tions, and the study of exotic quantum phases arising 
from the long-range anisotropic dipole-dipole interaction 
between polar molecules. 

Amongst the methods developed for the production of 
cold molecules (for an overview, see, for example, [IIHU] 
and references therein), those based on the deceleration 
of supersonic molecular beams are particularly suited to 
collision experiments pTMl7j : The final velocity of the 
decelerated particles is tunable over a wide range with a 
narrow velocity spread [HI [18] , the deceleration process 
is internal-state selective [111 HO] ■, and the initial phase- 
space density is preserved by operating the decelerator in 
a phase-stable manner [2TU26j . Depending on the char- 
acteristics of the species of interest, multistage Stark de- 
celeration [TB] , Rydberg-Stark deceleration [571 HH] , mul- 
tistage Zeeman deceleration [55H55] . or optical decelera- 
tion [34] may be best suited. 

Multistage Zeeman deceleration relies on the interac- 
tion between paramagnetic species and time-dependent 
inhomogeneous magnetic fields, and is therefore partic- 
ularly attractive for open-shell systems such as molec- 
ular radicals and metastable atoms and molecules. In 
this article, we present results on the deceleration of 
a supersonic beam of neon atoms in the metastable 
{lsf{2sf{2pf{isf state using a 91-stage Zeeman 
decelerator. Because of their favorable ratio of magnetic 



moment to mass and the straightforward detection of 
beams of metastable atoms, light metastable rare-gas 
atoms are well suited to multistage Zeeman decelera- 
tion, and beams of such atoms have been used by several 
groups to characterize their Zeeman decelerators [501^55] . 
Moreover, cold decelerated beams of metastable rare-gas 
atoms are of interest in the study of Penning ionization 
processes [35l|36| with a high energy resolution at low col- 
lision energies [371 13H] • Metastable neon has been laser 
cooled and trapped [321 ■ Multistage Zeeman decelera- 
tion also provides a means to prepare atomic or molec- 
ular samples in a selected magnetic sublevel, which may 
be used in collision experiments and complement exist- 
ing techniques to prepare aligned and oriented samples 
(see [inillT] and references therein). 

This article describes the extension of our modular 
high-repetition-rate Zeeman decelerator design, which 
has previously been employed for deceleration and trap- 
ping of hydrogen and deuterium atoms [551 H^Hii] , to 
the deceleration of heavier species, and also presents the 
characterization of the velocity distribution of the decel- 
erated samples. Characteristic features of this deceler- 
ator, next to its modular nature, are its high repetition 
rate (10 Hz) and the solenoid design which enables one to 
easily calculate magnetic field distributions. The latter 
advantage is of crucial importance for particle trajectory 
simulations and their use in the optimization of the de- 
celeration process. An introduction into design consider- 
ations relevant to multistage Zeeman deceleration will be 
given first, followed by a detailed description of the indi- 
vidual deceleration stages, their assembly into modules, 
and the integration of the modules into the multistage 
decelerator. The electronics circuits developed and used 
to pulse high currents through the deceleration solenoids 
with short rise and fall times are also described. 



II. EXPERIMENTAL SETUP 

In a multistage Zeeman decelerator, the longitudi- 
nal velocity of low-field-seeking particles in a supersonic 
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FIG. 1. (a) Schematic energy level diagram of the states of Ne relevant to the experiment, with wave numbers given relative 
to the (2p)^ ^So ground state [SU]. The metastable states and ^Po are indicated by thick red lines, (b) Zeeman effect in 
the state. The low-field-seeking state Mj = 2 with the largest magnetic moment is depicted in red. 



beam is reduced by converting kinetic energy into poten- 
tial energy through repeated application of pulsed inho- 
mogeneous magnetic fields 0. Typical laboratory mag- 
netic fields, which can be reliably produced and rapidly 
switched with conventional solenoids, are on the order of 
a few Tesla. These give rise to Zeeman shifts AWz/{hc) 
« 1 cm~^ for a magnetic moment of one Bohr magneton 
(/iB)- To give an example, the kinetic energy of a particle 
of mass number Nm propagating at a velocity of 400 m/s 
amounts to i?kin/(^c) ~ -^m x 6.7 cm~^. Many succes- 
sive deceleration stages are therefore needed to bring even 
small atomic or molecular systems to a standstill. More- 
over, a phase-stable operation of the decelerator is es- 
sential to maintain the initial phase-space density while 
decelerating parts of the supersonic beam to the desired 
final velocity [Hl[2i|2S]. 

A simplified energy level diagram of the metastable 
states of neon relevant to the experiment is shown in 
FIG. [l] The Zeeman shift of the metastable state is 
AWz = gjiiBMjB, with B the magnetic field and gj = 
1.50 |45j. At a magnetic field of 1.7 T, reached in our 
decelerator at a phase angle of 45° , the low-field-seeking 
state Mj = 2 used in the deceleration experiments is 
subject to a Zeeman shift of about 2.4 cm~^. 

A. Overview of the experimental setup 

A schematic representation of the experimental setup 
is shown in FIG. [2] (a). The main components are a 
pulsed source producing a supersonic beam of metastable 
neon atoms, the modular Zeeman decelerator, and the 
electronics driving the current through the solenoids to 
generate strong magnetic fields for deceleration. In the 
design of the decelerator and its electronics, the following 
aspects have been considered to make the experiment 
flexible and versatile: 

1. The decelerator was set up in a modular way to en- 
able the deceleration of different atomic or molec- 



ular systems. The number of deceleration stages 
required to fully decelerate a given species depends 
on its magnetic-moment-to-mass ratio. For exam- 
ple, 12 deceleration stages are needed for H 
24 deceleration stages for D [44], and around 90 
deceleration stages for metastable Ne. 

2. Collision and high-resolution spectroscopic experi- 
ments with pulsed supersonic beams are typically 
performed at repetition rates between 10 — 30 Hz. 
Ideally, the repetition rate of the decelerator should 
lie in this range, which places strong constraints on 
the cooling of the solenoids. 

3. To keep the number of deceleration stages low, the 
kinetic energy removed at each deceleration stage 
should be as large as possible. 

4. To avoid particle loss by nonadiabatic transitions at 
low field, a quantization field must be maintained 
along the entire trajectories of the decelerated par- 
ticles [T7|. 

5. The deceleration efficiency benefits from precise 
knowledge of the spatial and temporal shapes of 
the magnetic field pulses. In addition, "clean" 
magnetic-field pulses simplify the calculation of de- 
celeration pulse sequences. Metallic materials close 
to the solenoids should therefore be avoided to min- 
imize the influence of eddy currents. Moreover, 
high-permeability materials, which can be used for 
field-enhancement purposes [301 US] i result in a rest 
magnetization that is difficult to quantify. 

6. Whereas electronic devices for the fast switching of 
high voltages are commercially available, dedicated 
fast switching electronics for pulsing high currents 
through inductive loads are required. 

7. A high vacuum should be maintained in the decel- 
eration beam line to avoid losses by collisions with 
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FIG. 2. Schematic representation (not to scale) of the experimental setup. All distances are given in mm. (a) The modular 
decelerator with its main components: the metastable atom source, the evacuated volume with the tower platforms Tl-Tll, 
the towers, and the deceleration modules. The metastable atoms are detected after an adjustable free-propagation distance 
to the micro-channel-plate (MCP). (b) Dimensions of the deceleration solenoids and their assembly in a deceleration module, 
(c) Interface between the water-cooled deceleration modules and the pumping towers. The solenoids in the tower are not 
included for clarity. 



the background gas and with particles not accepted 
by the decelerator [S]. 

The detection of metastable atoms or molecules in a 
beam experiment is straightforward: If the internal en- 
ergy of the metastable state exceeds the work function 
of the micro-channel-plate (MCP) detector, an electron 
current is generated upon impact on the MCP. With an 
internal energy of 16.62 eV, the state of Ne clearly 
fulfills this requirement. Depending on the specific mea- 
surement, the MCPs used in the experiments described 
here were assembled in a standard chevron configuration 
for the measurement of time-of-flight (TOF) profiles, or 
additionally equipped with a phosphor screen for imag- 
ing the radial distribution of the decelerated atoms with a 
CCD camera. Compared to detection by photoionization 
with a pulsed laser, which was used in the Zeeman decel- 
eration of H and D [JT], this detection method has the 
advantage of giving a complete time-of-flight trace with- 
out the need of scanning the delay time of the laser with 
respect to the nozzle opening time. Hence, data acquisi- 
tion is fast, and a wide parameter range can be sampled 
in a short time. This advantage is particularly important 
during the initial optimization of the decelerator and the 
deceleration pulse sequences. 



B. Metastable- neon source 

Neon is expanded from a reservoir at a stagnation pres- 
sure of 5 bar into vacuum through the 250 |J.m diameter 
orifice of a pulsed Even-Lavie valve [?5^ . To increase the 
density in the beam and reduce the initial speed of the 
atoms, the body of the valve is thermally connected to a 
liquid-nitrogen reservoir by copper braids and held at a 
temperature of about 125 K. To populate the metastable 
•^P2 state, a discharge is operated in the expansion re- 
gion [49] . The distance between the nozzle orifice and 
the discharge is kept as short as possible to produce the 
metastable atoms in a high-density region, which was 
found to be efficient and also cools the beam in the ex- 
pansion. For this purpose, a thin metal tip, to which a 
DC voltage of -t-180 V is applied, is placed 1 mm away 
from the ceramic front plate of the valve body, with the 
grounded metal poppet of the valve acting as the second 
discharge electrode. The tip is radially displaced from 
the symmetry axis of the nozzle to avoid obstruction of 
the supersonic beam. To efficiently operate the discharge 
at this low voltage and minimize heating of the beam, the 
discharge is seeded with thermal electrons continuously 
emitted from a filament (Agar Scientific A054C, 0.2 mm 
tungsten wire, AC current up to 6 A rms) placed down- 
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stream from the nozzle. 

To characterize the supersonic beam, we have recorded 
the time-of-fiight distribution of the metastable neon 
atoms impinging on the MCP at the exit of the decel- 
erator, about 1.8 m from the nozzle. From such time-of- 
flight measurements, the mean velocity of the metastable 
neon atoms was determined to be 580 — 590 m/s. Fur- 
ther properties of the discharge source can be determined 
by threshold photoionization spectroscopy using a pulsed 
dye laser. After propagation through the decelerator, 
only the metastable and ^Pq states are expected 
to be populated in the beam. Scanning the laser over 
their ionization thresholds at 39111.11 cm~^ (^-Fo) and 
39887.91 cm-i (^Ps) [5D], and recording the ion yield, 
we found a ratio of population in the state to that 
in the "^Pq state of 5:1, which corresponds to the degen- 
eracy factors of these energy levels. Because the state 
'^P2, Mj — 2 has the largest magnetic moment it is the 
most convenient state for deceleration, and the deceler- 
ation pulse sequences were optimized for its properties. 
Assumining an equal initial population of all magnetic 
sublevels, the atoms in the Mj = 2 state represent only 
one sixth of the population of metastable atoms. In the 
measurements, neon gas with natural isotopic composi- 
tion was used, the most abundant isotopes being ^°Ne 
(90.5%) and ^^Ne (9.3%). In most experiments, the de- 
celeration pulse sequences were generated for ^°Ne. The 
occurrence of two different isotopes in the beam, how- 
ever, offers the possibility of investigating the selectiv- 
ity of multistage Zeeman deceleration on the magnetic- 
moment-to-mass ratio. This was done by applying decel- 
eration pulse sequences generated for ^°Ne and ^^Ne, and 
detecting the cold decelerated atoms in a mass-selective 
manner. 



C. Multistage Zeeman decelerator 

The deceleration solenoids are at the heart of the mul- 
tistage Zeeman decelerator. In their design and optimiza- 
tion, the following factors have been considered. 

1. The open area of the decelerator, determined by 
the free inner diameter of the solenoids, should be 
matched to the profile of the supersonic beam. The 
dimensions are therefore similar to those of the elec- 
trodes in a multistage Stark decelerator. 

2. The solenoids should produce a large decelerating 
field gradient on axis and provide transverse focus- 
ing of the beam. 

3. The length of the solenoids should be kept as short 
as possible to maximise the energy removed per 
unit length. 

4. The length of the solenoids should be matched to 
the initial speed of the supersonic beam and to the 




FIG. 3. The axially symmetric magnetic field produced by 
the 7.2 mm-long deceleration solenoid, centered at z=0 and 
operated at a current of 300 A. On the symmetry axis of the 
solenoid, the barrier has a height of 2.2 T. The thick curves 
on the surface at z = mm, z = —5 mm, and z = — 10 mm 
highlight the concave shape of the field near the center of 
the solenoid, which gives rise to radial focusing forces, and 
the convex shape at larger distances from the center, where 
forces become radially defocusing. The longitudinal extension 
of the solenoid is indicated by the thick bar on the z-axis. 



achievable switching times of the magnetic decel- 
eration fields. For efficient phase-stable decelera- 
tion, the beam should not penetrate too far into 
the solenoids during the switch-off time of the field. 

The individual deceleration solenoids were made of 
enameled wire having a copper-core diameter of 400 |^m 
(Elektrisola XHTW, ABP15, grade IB). The wire was 
wound in four layers with 62 windings in total, with a 
free inner solenoid diameter of 7 mm, an outer diame- 
ter of 10.4 mm, and a length of 7.2 mm. The length 
of these solenoids was chosen as a compromise between 
achieving the maximum on-axis field for a given number 
of windings per unit length at this inner diameter while 
at the same time maximizing the deceleration that can 
be achieved per unit length. Electrically, the solenoids 
can be characterized by an inductance of 21 |j.H and an 
ohmic resistance of about 0.24 O at room temperature. 
To improve the mechanical stability of the solenoids, a 
supporting structure is produced by curing the enamel 
coating. 

The magnetic field distribution of a solenoid calculated 
for a typical current of 300 A, the maximum pulsed cur- 
rent rating of the IGBT (insulated-gate bipolar transis- 
tor) switches used, is depicted in FIG. [3] The solenoids 
create a magnetic field barrier of 2.2 T along the su- 
personic beam propagation axis, which is used for de- 
celeration. Near the center of the solenoid, the field is 
concave, which gives rise to radial focusing forces. At 
larger distances from the center, the field is convex, and 
forces in the radial direction are defocusing. This alter- 
nation plays an important role in the transverse dynam- 
ics in the multistage Zeeman decelerator as discussed by 
Wiederkehr et al. |26) . 
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When a current of 300 A is pulsed through a solenoid, 
an electrical power of 20 kW is dissipated. At a rep- 
etition rate of the decelerator of 10 Hz and a typical 
pulse duration of 50 |i.s, a heat load of 10 W per solenoid 
must be removed by the cooling system. To efficiently 
cool the solenoids and simplify the construction, 12 in- 
dividual solenoids are assembled into one module as 
shown in FIG. [2] Within such a module, the super- 
sonic beam propagates inside a quartz glass tube (in- 
ner diameter 5 mm, wall thickness 1 mm) representing 
the evacuated beamline. The solenoids, separated by ce- 
ramic (BNP2) spacers at a center-to-center distance of 
10.7 mm, are glued on the outside of this glass tube with 
a heat-conducting epoxy (Aremco-Bond 860) . The epoxy 
prevents mechanical distortions when high currents are 
pulsed and protects the outermost layer of the solenoid 
from abrasion by the circulating cooling water. The glass 
tube assembled with the solenoids is then inserted into 
a plastic casing which also houses electrical feedthroughs 
to supply the current to the solenoids. This design al- 
lows the immersion of the solenoids in circulating cool- 
ing water to efficiently remove the heat, such that the 
decelerator can be operated at a high repetition rate of 
10 Hz. 

Before the supersonic beam enters the decelerator, it 
passes a 2 mm diameter skimmer placed 10 cm down- 
stream from the nozzle. The nozzle is positioned in the 
source vacuum chamber with good optical access, which 
is necessary when the species to be decelerated are pro- 
duced by either photoexcitation or photolysis from a pre- 
cursor molecule. The size of the skimmer orifice is chosen 
so as to match the transverse profile of the supersonic 
beam to the radial acceptance of the decelerator. To re- 
duce the gas load from the source in the decelerator and 
avoid losses by collisions with the background gas |8] , the 
beam first propagates through a tower extending above 
an evacuated volume (see FIG. [2]). The large vacuum 
chamber, covering the full length of the decelerator, is 
evacuated by two 500 1/s turbo molecular pumps to a 
pressure in the low lO"'' mbar range. 

For the deceleration of Ne ("^^2)7 five modules of 
144 mm length, each equipped with 12 solenoids, and 
one additional 348 mm-long module with 31 deceleration 
solenoids were used. As shown in FIG. [2j the mod- 
ules are attached with viton 0-ring seals to the towers. 
Therefore, modules can be added or removed easily to 
adapt the length of the decelerator to different species. 
The 12-solenoid modules are similar to the ones used for 
the deceleration and trapping of H and D atoms, and 
were designed so as to permit deceleration of H atoms to 
a standstill '5i|lHllj. 

Interrupting the sequence of deceleration stages by the 
60 mm-long propagation distances through the towers is 



When the 348 mm-long module is used, one tower is skipped, 
and the length of the module is given by twice the tower spacing 
plus the tower length. 



necessary to maintain a high vacuum in the decelerator, 
but it introduces additional technical difficulties. In the 
case of neon, only atoms in the ^P2,Mj = 2 state are 
efficiently decelerated. To avoid nonadiabatic transitions 
of the low-field-seeking states near zero field, which would 
cause particle loss O |47] , it is important to maintain a 
quantization field along the complete trajectory through 
the decelerator. Within the modules, this is achieved by 
switching on the field of the next solenoid just before the 
field of the active solenoid is switched off. In each tower, 
two larger solenoids (free inner solenoid radius 6 mm, six 
layers with 11 windings each, 1 mm diameter enameled 
copper tubing with a wall thickness of 0.2 mm, center- 
to-center spacing of the solenoids 28 mm) , depicted in 
FIG.[2j define the quantization axis H^. The adjustment 
of the currents of these solenoids can also be used to 
manipulate the transverse motion of the atoms. 



D. Electronics driving the deceleration solenoids 

To operate the multistage Zeeman decelerator in an 
efficient way, the particles must be prevented from mov- 
ing over a long distance during the switch-off time of the 
field ^26j. Therefore, the decelerating magnetic field of 
each stage must be switched off as rapidly and completely 
as possible. 

In FIG. |4j we show a diagram describing the opera- 
tion principle of the electronic circuit used to achieve the 
sub- 10 |as rise and fall times of the current pulses ap- 
plied to the deceleration solenoids. The circuit relies on 
a high-voltage supply for fast switching together with a 
low- voltage supply to maintain a constant current during 
the hold phase. To control the currents, insulated-gate 
bipolar transistor (IGBT) switches with integrated free- 
wheeling diodes (Semikron, SKM200GB176D) are used, 
which permit operation with pulsed currents of up to 
300 A. 

To rapidly ramp up the current through the solenoid, 
the IGBT switches K and S are closed, and the voltage 
Uk is applied to the solenoid. The current through the 
solenoid increases as Igoi = lo (1 — exp(— i/r)) « lot/r 
(Iq = Uk/R, t = L/R, L = 21 |xH, R = 0.24 n), the 
linear approximation being valid at short times. The 
current through the solenoid is monitored by the voltage 
drop across the shunt resistor Rg and fed to a compara- 
tor. When a preset value of the current has been reached 
at the end of the switch-on phase (1) (blue color code in 
FIG. [4]), switch K is opened. The current commutes to 
branch (2) (green color code in FIG. |4]) and is held con- 
stant because of the high energy stored in the capacitor 
Ch. 

To switch off the magnetic field produced by the 
solenoid, switch S is opened at the end of the hold phase 
(2). The decay of the magnetic field in the solenoid is ac- 
companied by the build-up of an induction voltage which 
rises to a potential Uk+Ud (Ud: forward voltage drop of 
diode Dd) and opens a path (3) (red color code in FIG.|4]) 
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FIG. 4. (a) Operation principle of the electronic circuit used 
to generate pulsed magnetic fields with short rise and fall 
times. The deceleration solenoid is represented by its equiva- 
lent resistance and inductance inside the dashed frame. Uk: 
kick voltage supply (800 V); Ck: storage capacitor (50 viF); 
Uh: hold voltage supply (100 — 120 V); Ch: storage capac- 
itor (26.4 mF); Dh, Dd: hold diode and decay diode; S, K: 
IGBT (insulated-gate bipolar transistor) switches; Rs: shunt 
resistor (10 mil). The current flow during the switch-on (1), 
hold (2), and switch-off (3) phases is indicated by the blue, 
green, and red arrows next to the wires. The nominal values 
of Ck and Ch are chosen to supply the current pulses to six 
deceleration solenoids (see FIG. [6|. (b) Schematic temporal 
profile of the current through the solenoid, (c) Logical states 
of the IGBT switches during the three switching phases (0: 
open, no current flow; 1: closed). 
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FIG. 5. Magnetic field at the center of a deceleration solenoid. 
The field was obtained by integrating the voltage induced 
across a pickup coil (1.4 mm free inner diameter, 15 wind- 
ings of a 0.1 mm-diameter enameled copper wire, L=250 nH) 
placed in the bore of the solenoid. 



for the solenoid current into the storage capacitor Ck. 
The current during the switch-off phase decreases with a 
time-reversed temporal behavior compared to the switch- 
on phase, and stops when all of the energy stored in the 
solenoid has been either transferred into capacitor Ck or 
dissipated as heat in the circuit. 

The temporal magnetic-field profile of a deceleration 
solenoid as measured with a pickup coil is shown in 
FIG. [5| At a kick voltage of Uk=800 V, a rise and 
fall time of the 2.2 T magnetic field of 8.5 |i.s is deduced 
from this measurement. For typical durations of the hold 
phase of around 30 |j.s, the magnetic field strength re- 
mains constant, with a negligible decrease (< 2%) re- 
sulting from a slow discharge of the hold capacitor Ch. 
The undershoot at the end of the switch-ofT phase is at- 
tributed to the properties of diode Dd and has to be com- 
pensated to maintain a well-defined quantization field 
throughout the decelerator. This is achieved by apply- 
ing a current to the adjacent deceleration solenoid. The 
weak oscillations following the undershoot, however, are 
an artifact caused by capacitive coupling to the pick-up 
coil. 

The switching electronics are assembled into units 
capable of driving one deceleration module with 12 
solenoids. Each of these electronics units consists of 
two drivers labeled A and B. Both drivers have their 
own, independent power electronics schematically drawn 
in FIG. [6] (a) and a control and interface circuitry as 
shown in FIG. [6] (b). This subdivision into two drivers 
A and B results from the necessity to maintain a quan- 
tization field at the position of the decelerated particles 
at all times, which is achieved by temporally overlapping 
the current pulses applied to adjacent solenoids as illus- 
trated in FIG. El To minimize the total number of IGBT 
switches, a single IGBT switch K is used to pulse current 
through six of the deceleration solenoids in one module. 

The individual electronics units are connected to 
a general-purpose pattern generator (Spincore Pulse- 
Blaster PB24-100-512) via a bus system, allowing for fiex- 
ibility in the generation of deceleration pulse sequences. 
One electronics unit contains two identical interface and 
control circuits, each of them servicing one of the two 
drivers A and B. The bus interface decodes the signals 
from the pattern generator, and extracts the addresses 
of the electronics unit and the solenoid to be energized. 
To control the temporal current profile applied to the 
solenoid, the voltage Us across the shunt resistor Rs is 
monitored, and IGBT switch K is opened if the preset 
value of the current (usually 300 A) is reached. 

In addition to controlling and switching the power elec- 
tronics, the control circuitry provides protection and in- 
terlock functions. The protection circuit inhibits gate 
drive to the IGBT if the repetition rate is too high, if 
current pulses longer than a specified maximal duration 
are applied, or if the interval between successive pulse 
commands is too short. In addition, the status output of 
the integrated monitor electronics of the IGBT driver is 
processed, thus inhibiting gate drive to the IGBT in case 
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FIG. 6. (a) Schematic diagram of the power electronics driver A used for pulsing the current through deceleration solenoids 
A1-A6 in a deceleration module consisting of 12 solenoids. Uk(h)5 Ck(h): kick (hold) voltage supply and storage capacitor as 
defined in FIG. |4] K: control input for kick voltage switch; S1-S6: control input for the switches used to energize one of the 
deceleration solenoids A1-A6 at a time. The switch action corresponds to that of S in FIG.|4] (b) Schematic diagram of the 
interface and control circuitry for driver A. The circuitry bounded by the dashed line is realized as a synchronous design with 
CMOS programmable logic devices (CPLD, Lattice M4A5), resulting in a time resolution of 100 ns. 



of an overcurrent condition. For monitoring purposes, 
the state of the electronics unit is displayed using sta- 
tus indicator LEDs, and the solenoid current measured 
across the shunt resistors is supplied to an analog inter- 
face. 

To pulse current through the tower solenoids, a simi- 
lar power electronics and control circuitry is used. How- 
ever, the higher inductance (50 |i.H per solenoid) of these 
solenoids imposes longer switching times and necessitates 
the use of larger storage capacitors to provide a suffi- 
ciently high energy reservoir. Therefore, a single-channel 



electronics unit supplies current to the two solenoids of 
each tower which are connected in series. 



III. RESULTS AND DISCUSSION 

The deceleration process was characterized by a combi- 
nation of time-of-flight and imaging techniques. To gen- 
erate a deceleration pulse sequence, the synchronous par- 
ticle with a given initial velocity was propagated along 
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FIG. 7. (a) Alternating connection of the deceleration 
solenoids in one deceleration module to the power electron- 
ics drivers A and B. (b) To maintain a quantization field at 
all times during deceleration, the current through solenoid 
N+l is switched on before switching off the current through 
solenoid N. 



the axis of the decelerator in a one-dimensional trajec- 
tory simulation taking into account the finite rise and fall 
times of the decelerating fields. As explained in we 
define the phase angle ipo in a deceleration sequence as 
the reduced position of the synchronous particle at the 
time the switch-off of the field is initiated. When a mea- 
surement is referred to as having been performed at a 
phase angle (j)o, this means that all deceleration stages 
throughout the decelerator were operated with the same 
phase angle (j}Q. Parameters which were systematically 
varied in the measurements include the phase angle and 
the initial velocity of the synchronous particle. 



A. Time-of-flight distribution of decelerated atoms 

Time-of-flight distributions recorded for different set- 
tings of the decelerator are compared in FIG. |8]to a time- 
of-flight distribution obtained with the decelerator inac- 
tive (gray trace). In these measurements, time zero corre- 
sponds to the trigger of the pulsed valve. For a supersonic 
beam at high Mach number, the velocity distribution can 
be approximated by a Gaussian distribution. 

When the decelerator is turned on, the appearance of 
the main time-of-flight peak changes qualitatively (see 
lowest black trace in FIG. [8] recorded at a phase angle 
of 0°). The rich substructure of this peak is caused by 
transverse focusing and guiding of the metastable atoms 
by the pulsed magnetic fields [HI S?]- In this time-of- 
flight profile, an additional peak is observed at a later 
time (~ 4.6 ms), corresponding to neon atoms that are 
decelerated in a phase-stable manner at (f)o =0°. It is, 
however, not possible to directly extract the velocity of 
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FIG. 8. Time-of-flight distribution of metastable neon atoms 
recorded with the MCP detector placed 370 mm downstream 
of the decelerator exit. The gray trace was measured with 
the multistage decelerator inactive. For the other traces, the 
deceleration sequence was calculated for an initial velocity of 
the synchronous particle of 580 m/s, and the phase angle was 
varied between 0°-55° as indicated next to the trace. For 
clarity, the traces recorded with phase angle f Q°-55° are ver- 
tically offset, and only a time window centered at the arrival 
time of the decelerated atoms is shown. 



the decelerated atoms from their arrival time with re- 
spect to the time of flight of the undecelerated atoms, 
because the change in energy and velocity per time inter- 
val is not constant during the propagation of the atoms 
through the decelerator. Nevertheless, when all decel- 
eration stages are operated at the same phase angle, a 
later arrival time of the decelerated packet corresponds 
to a lower flnal velocity. Comparing the signal of the 
decelerated-atom peaks in FIG. [8] to the signal of atoms 
with the decelerator inactive, it should be noted that all 
five magnetic sublevels of the state as well as the '^Pq 
state contribute to the latter, whereas only atoms in the 
state ■^P2, Mj — 2 contribute to the peak of decelerated 
atoms. 

When a deceleration sequence with an increased phase 
angle is applied, a larger amount of kinetic energy is re- 
moved per deceleration stage, and the particles within the 
phase-stable region around the synchronous particle are 
decelerated to lower velocities. This behavior can be ob- 
served in the time-of-fiight distributions shown in FIG.|8) 
where the arrival time of the decelerated atoms increases 
with increasing phase angle. At the same time, a one- 
dimensional model of a multistage Zeeman decelerator 
predicts a reduction of the phase-stable volume accepted 
by the decelerator when the phase angle is increased ; 
one might thus expect the largest signal of decelerated 
atoms at 0° phase angle. The intensity of the decelerated- 
atom peak shown in FIG. [9) however, does not show such 
a monotonic behavior. Instead, we find that the decel- 
erator appears to operate most efficiently at a phase an- 
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FIG. 9. Intensity of the decelerated atom signal deter- 
mined from the area of the decelerated-atom peaks in the 
time-of-flight distributions. The deceleration sequences were 
generated for an initial velocity of the synchronous particle 
of 580 m/s. 



gle of 40°. A careful analysis of the three-dimensional 
particle dynamics in the multistage Zeeman decelerator 
reveals strong transverse effects not accounted for in the 
one-dimensional model. While transverse focusing and 
defocusing effects reduce the phase-stable volume at low 
phase angle, the finite rise and fall times of the pulsed 
magnetic fields in the deceleration solenoids tend to com- 
pensate this to some degree and increase the phase-stable 
volume. Altogether, three-dimensional trajectory simu- 
lations demonstrate that with our solenoid configuration 
the multistage Zeeman decelerator has the largest ac- 
cepted phase-stable volume («10mm'^x200 (m/s)"^) for 
phase angles in the range 30°-45° [26], in agreement with 
our experimental findings (see Sections III B and III C I . 
The present analysis further confirms that the results 
of the phase-stability analysis performed for D (1^5*1/2) 
atoms can be transferred to the deceleration of heavier 
species in significantly longer decelerators. 

At the base pressure of the decelerator of 
2 X 10^^ mbar, losses by collisions with the back- 
ground gas are negligible. Experiments carried out at 
higher pressure by introducing nitrogen into the vacuum 
chamber through a needle valve indicate that colli- 
sional losses become significant at pressures exceeding 
2 X 10"^ mbar. 



FIG. 10. Time-of-flight distributions of metastable neon 
atoms recorded with 200 mm, 370 mm, and 505 mm of free 
propagation from the exit of the decelerator to the MCP de- 
tector. The deceleration pulse sequence was generated for an 
initial velocity of 580 m/s and a phase angle of 40°. For clar- 
ity, only a portion of the time-of-flight distribution centered 
at the arrival time of the decelerated atoms is shown. The 
solid curve is a fit of S{Li,t) dt, with S{L, t) dt defined by 
Eq. ||8|, to these data. 



distribution by free propagation. By recording time-of- 
flight distributions for several distances of free propaga- 
tion (200 mm, 370 mm, and 505 mm) between the exit of 
the decelerator and the MCP detector, it is possible to 
reconstruct the velocity distribution. 

The timc-of-flight profiles of decelerated neon atoms 
recorded at a phase angle of 40° for these three differ- 
ent free-propagation distances are compared in FIG. [TO] 
Similar profiles were recorded for several phase angles 
beween 0° and 55°. The increasing temporal spread of 
the packet of decelerated atoms during its propagation to 
the detector is caused by velocity dispersion. To extract 
the longitudinal temperature from the time-of-flight pro- 
files, two methods were used. The first method, based on 
Gaussian fits to the observed time-of-flight distributions, 
offers the advantage of being robust and simple, and en- 
ables rapid diagnostics of the deceleration process during 
data acquisition. The second, more accurate method, 
relies on a complete model of the particle propagation 
including spatial and velocity distributions. For the anal- 
ysis, all traces were individually normalized to the signal 
of decelerated atoms to account for small variations in 
the operation conditions of the metastable neon source. 



Determination of the longitudinal velocity 
distribution 



The direct detection of metastable neon atoms im- 
pinging on the MCP offers a simple way to experimen- 
tally characterize the velocity distribution of the beam: 
the velocity distribution is mapped onto an arrival-time 



1. Velocity of decelerated atoms 
For a preliminary analysis, a fit of Gaussian functions 
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FIG. 11. Measured final velocities of decelerated neon atoms 
determined from their arrival time at the detector with dif- 
ferent distances of free propagation. In the case of the lowest 
final velocity, 106 m/s, more than 95 % of the initial kinetic 
energy of the beam was removed in the 91-stage decelerator. 



was made to the individual peaks of decelerated atoms 
in the time-of-flight signal recorded for different free- 
propagation distances. The mean velocity v and the ve- 
locity spread CTi, were then extracted from an analysis 
of the peak position tc{L) and the width <Tt{L) of the 
decelerated-atom peak in the time-of-flight distribution 
for each free-propagation distance L. The mean velocity 
of the decelerated atoms was determined by fitting 



tc{L) 



V 



(2) 



to the data, where is the time at which the free prop- 
agation starts (i.e., the time the decelerated atoms have 
reached the exit of the decelerator). The velocity spread 
of the decelerated atoms was then obtained by fitting 



2. Modeling the time-of-flight distribution 

To model the arrival-time distribution at the detec- 
tor, the longitudinal velocity distribution of the packet 
of atoms decelerated to a velocity v is described by a 
Gaussian distribution 



f{v)dv 



1 



2tT(7, 



exp 



{v — v)' 
2al 



du, 



(4) 



where cr^ is the variance in the final longitudinal velocity. 
Free propagation over a distance L maps the velocity 
distribution /(w) dw to an arrival-time distribution at the 
detector. 



S{L,t) dt = f{Lt-^)Lt-^dt. 



(5) 



The longitudinal extent of the cloud of decelerated atoms 
at the beginning of the free-propagation phase can be 
approximated by a Gaussian distribution. 



g{z) dz 
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The arrival-time distribution is obtained by the convolu- 
tion 



S{L,t) dt 



/( 



t 



')Lt-^g{z) 



dz dt, 



(7) 



which, in the case of two Gaussian distributions, yields 

{L-vtf 



S{L,t) dt = 
1 Ltal 



2^{t^ a2 + a2)3/2 



exp 



dt. (8) 



This expression has the same form as Eq. ([T]), which jus- 
tifies the use of Eq. ([2| and Eq. ([3| in the first analy- 
sis, but links the observed time-of-flight profiles to both 
the spatial and velocity distributions of the decelerated 
atoms. 



(3) 



3. Velocity distributions from time-of-flight measurements 



to the data, where CTz is the width of the Gaussian packet 
at the exit of the decelerator. In FIG. 11 the velocities of 
the decelerated atoms for different phase angles are com- 
pared. These velocities were obtained by analysing mea- 
suremets for three different free-propagation distances 
and correspond to the final velocities for which the decel- 
eration pulse sequences were optimized to within 1 m/s. 
This indicates that the multistage decelerator and the de- 
celeration process are very well characterized and atom 
bunches of a desired final velocity can be reliably gen- 
erated. The analysis also confirms that a higher phase 
angle results in a lower velocity of the decelerated atoms, 
as expected from simple considerations. 



To reconstruct the velocity distribution of the de- 
celerated atoms, ^jS'(Li,t)dt {i = 1,2,3) was fitted 
to the time-of-fiight distributions measured with free- 
propagation distances Li = 200 mm, L2 = 370 mm, and 
L3 = 505 mm as shown in FIG. [TOj This procedure is 
more sensitive to changes of the peak shape with propa- 
gation than fitting individual Gaussian functions to the 
peaks in the time-of-fiight traces. Thus, a good agree- 
ment between the fit and the experimental data, as shown 



in FIG. 10 is only obtained if the model functions used 
for the velocity distribution and the initial spatial distri- 
bution are good approximations to the real phase-space 
distribution of the decelerated atoms. The parameters 
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FIG. 12. (a) Velocity spread of the decelerated atoms determined by fitting the model function S{Li,t) dt to the time- 
of-flight distributions. The data points in a group of measurements at the same phase angle but different initial velocities 
are slightly displaced horizontally from each other for clarity. The error bars estimate the accuracy of the fit in reproducing 
the data. The temperature is determined from the velocity spread using the relation (b) Velocity spread and 

(c) phase-space distribution of the decelerated atoms obtained from three-dimensional numerical particle-trajectory simulations. 



which can be obtained from the fit are the velocity of 
the decelerated atoms f dccoi and their velocity spread a^, . 
Although the model function Eq. ([s]) depends on the lon- 
gitudinal position spread at the exit of the decelerator, 
the fits arc rather insensitive to this parameter because 
the temporal spread of the packet in the time-of-flight 
distribution is dominated by velocity dispersion, even for 
the shortest distance of 200 mm. Therefore the longitu- 
dinal spread at the exit of the decelerator obtained from 
trajectory simulations, which was found to decrease from 
az — 2.5 mm at 0o = 0° to cr^ = 1 mm at 0o — 50°, was 
used as a fixed input parameter to the fit. 

In FIG. |12[ the experimentally determined depen- 
dence of the velocity spread of the decelerated atoms 
on the phase angle is compared to the outcome of 
three-dimensional numerical particle-trajectory simula- 
tions. The good overall agreement indicates that the sim- 
ulations reproduce the main features of the deceleration 
process adequately. For phase angles in the range of 30°- 
50°, where the decelerator is operated most efhciently 
(see FIG. [o] and Wiederkehr et al. [Mj), the longitudi- 
nal temperature T of the decelerated atoms determined 



from ksT = ^mcr^ is found to be only 10 mK. This 
is less than typical depths of magnetic traps of around 
100 mK [nillallll], so that efficient trap loading follow- 
ing the procedure described in [331 [33] can be envisaged. 

The temperature of the decelerated atoms is found to 
decrease with increasing phase angle. Following a sim- 
ple model of the Zeeman decelerator, which assumes in- 
finitely short rise and fall times of the decelerating fields, 
one expects the final temperature to directly refiect the 
phase-space acceptance of the decelerator. Hence, at a 
fixed phase angle, it should be independent of the initial 
velocity of the atoms. The simulations, however, indi- 
cate an influence of the initial velocity of the synchronous 
particle, for which the deceleration pulse sequence is cal- 
culated, on the temperature of the decelerated atoms. 
This effect is caused by the finite rise and fall times of 
the decelerating magnetic field. Atoms which move at a 
higher velocity propagate further into the solenoid dur- 
ing the switch-off time of the field. Thus, they experience 
a larger effective phase angle, which reduces the phase- 
space acceptance of the decelerator, and results in a lower 
temperature. 
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For comparison to previous work [26J , the longitudinal 
phase-space distribution of decelerated atoms obtained 
from trajectory simulations and the separatrix obtained 
from a one-dimensional phase-stability analysis are pre- 
sented in FIG. 12 (c). The phase-stable area in this dia- 
gram is found to be smaller than the one-dimensional sep- 
aratrix at each phase angle plotted. At low phase angles 
(here shown for 0°), a reduction of the number of par- 
ticles near the phase-space position of the synchronous 
particle, i.e., in the central region of the separatrix, is 
observed. As explained in [55], these differences between 
the predictions of the one-dimensional model and the 
phase-stable volume of a "real" multistage Zeeman de- 
celerator are caused by an interplay between the three- 
dimensional magnetic field distribution and the finite rise 
and fall time of the field. At low phase angles, the con- 
vex, defocussing nature of the magnetic-field distribution 
at large distances from the center of the solenoids cannot 
be compensated any more and a loss of particles results. 
The interruption of the otherwise regular array of decel- 
eration solenoids by the towers can cause an additional 
reduction of the phase-stable volume. AUtogether, the 
analysis demonstrates that the results of the extensive 
study of phase-stability performed for D atoms can 
be transferred to the description of a modular multistage 
Zeeman decelerator for other species. 



C. Imaging decelerated neon atoms 

The determination of the longitudinal velocity distri- 
bution of decelerated atoms by the method discussed in 
Section [ill B| was greatly faciliated by the direct detec- 
tion of the metastable neon atoms impinging on the ac- 
tive surface of the MCP detector. The timc-of-flight mea- 
surements of the longitudinal velocity distribution can be 
extended in a straightforward way to also determine the 
transverse spatial and velocity distributions of the decel- 
erated atoms. For this purpose, a MCP equipped with 
a phosphor screen behind the active layers was installed, 
and the phosphor screen was imaged from behind with a 
CCD camera. 

In FIG.[T3[ we show background-subtracted false-color 
images of decelerated metastable neon atoms obtained 
by this method. Such images were recorded for decelera- 
tion pulse sequences with a phase angle between 0° and 
50° , and for free-propagation distances of 200 mm and 
400 mm from the exit of the decelerator to the MCP. 
These images serve the purpose of characterizing the 
beam temperature and aid in the alignment of the de- 
celerator. 

In the images shown in FIG. [T3[ a displacement of the 
decelerated atom cloud from the reference position, given 
by the point of intersection of the decelerator axis and 
the MCP detector surface, is apparent. The shape of 
the atom cloud is also slightly distorted, the effect be- 
ing most pronounced at small phase angles. Both effects 
were found to sensitively depend on the transverse align- 
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FIG. 13. Background-subtracted false-color images of 
metastable neon atoms decelerated to final velocities of 
220 m/s and 180 m/s at a phase angle of 30° and 40°, re- 
spectively. L is the distance of free propagation between the 
exit of the decelerator and the MCP. The camera exposure 
time of 100 — 200 ]is was matched to the arrival-time spread 
of the decelerated atoms. The black cross indicates the point 
of intersection of the decelerator axis and the MCP detector 
surface. 



ment and tilting of the deceleration modules with respect 
to each other, and the images shown in FIG. [T3| repre- 
sent the best compromise achieved when optimizing the 
alignment for the most relevant phase-angle range of be- 
tween 30° and 40°. The asymmetry and displacement is 
likely caused by the limited precision in the positioning 
of the current-carrying solenoids. Since these effects be- 
come visible only at long distances from the exit of the 
multistage Zeeman decelerator, we do not expect them 
to hamper the use of the decelerator in high-resolution 
spectroscopy, in collision experiments, or in trapping ex- 
periments, which would ideally be performed at shorter 
distances. 

To obtain the transverse temperature and the trans- 
verse position spread of the packet of decelerated atoms 
at the exit of the decelerator, the transverse extension of 
the slow atom cloud was studied as a function of prop- 
agation distance. The analysis was performed for phase 
angles in the range 30°-45° in which the decelerator is 
operated most efficiently. Assuming a Gaussian distribu- 
tion, a transverse temperature of 5 — 10 mK can be de- 
termined in combination with a transverse spread at the 
exit of the decelerator of about 1 mm. These values are 
in agreement with the predictions of three-dimensional 
particle-trajectory simulations. This agreement demon- 
strates that the phase-space distribution of the deceler- 
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FIG. 14. Time-of-flight distribution of metastable neon 
atoms. The lower traces (A1-A3) were recorded with deceler- 
ation pulse sequences generated for ^''Ne, whereas the upper 
traces (B1-B3) were recorded with pulse sequences for ^^Ne. 
The initial velocity of the synchronous particle was 580 m/s 
(Al, Bl), 570 m/s (A2, B2), and 560 m/s (A3, B3), and a con- 
stant phase angle of 40° was used. The final velocities of the 
decelerated atoms are 191 m/s (Al), 161 m/s (A2), 129 m/s 
(A3), 248 m/s (Bl), 226 m/s (B2), and 202 m/s (B3). For 
clarity, the traces recorded with different deceleration pulse 
sequences are vertically offset. 



ated atoms can be completely reconstructed from a com- 
bination of timc-of- flight and imaging measurements. 



D. Isotope-selective multistage Zeeman 
deceleration of metastable neon 



The significant natural abundance of two stable iso- 
topes (2°Ne: 90.5%, ^^Ne: 9.3%) of neon opens up the 
possibility to investigate whether multistage Zeeman de- 
celeration can be used to separate species of different 
magnetic-moment-to-mass ratio. For this purpose, pulse 
sequences generated for a synchronous particle of each 
of the two isotopes were used, and the time-of-flight 
traces of decelerated metastable neon atoms were ana- 
lyzed. Typical results are displayed in FIG.[T4j Because 
the same amount of kinetic energy is removed from the 
synchronous particle per deceleration stage if the decel- 
erator is operated at constant phase angle, a larger veloc- 
ity change is expected for the lighter isotope. Comparing 
the time-of-flight distributions recorded for the same ini- 
tial velocity but different mass of the synchronous par- 
ticle (i.e., trace AI-A3 with B1-B3), the later arrival 
time of the decelerated atom peaks with the pulse se- 
quences generated for ^"Ne indicates a lower final ve- 
locity. This is confirmed by measurements with different 
free-propagation distances (see caption of FIG. 14 for the 
final velocities of the decelerated neon atoms determined 
experimentally) . 



However, these time-of-flight measurements alone do 
not allow one to disentangle the isotopic composition of 
the decelerated packet. To achieve this, a complemen- 
tary measurement based on photoionization was used, 
which allowed isotopic separation of the resulting pho- 
toions in a time-of-flight mass spectrometer. The results 
are displayed in FIG. [15] and show the time-of-flight dis- 
tribution of ions produced by pulsed-field ionization after 
photoexcitation of neon atoms from the state to Ry- 
dberg states of principal quantum number n = 24. The 
pulsed field also served the purpose of extracting the Ne"*" 
ions out of the metastable atom sample towards the MCP 
detector. 



In a first measurement shown in FIG. 15 (a), the decel- 
erator was inactive, and the delay of the laser pulse with 
respect to the valve opening time was set so as to pho- 
toexcite atoms at the maximum of the gas pulse (3.2 ms 
in FIG. 14 1. In this measurement, the ratio of the two 
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isotopes, as determined from the area of the correspond- 
ing peaks in the time-of-flight mass-spectrometer signal, 
was found to be in good agreement with that expected 
from their natural abundance. 

In a second measurement, a deceleration pulse se- 
quence generated for ^°Ne was applied, and the laser 
delay was set to detect atoms at the center of the de- 
celerated atom peak (4.95 ms in trace Al in FIG. 
the ion time-of-flight distribution shown in FIG. 
only ^°Ne is observed, and the signal at the expected 
arrival time of ^^Ne is reduced to the background level. 

In the last measurement, a deceleration pulse sequence 
generated for ^^Ne was applied. In this case, differences 
in the ion time-of-flight signal are observed depending on 
the time of photoexcitation within the decelerated atom 
peak in trace Bl in FIG. [W] In the ion time-of-flight 
signal shown in FIG. 15 (c), which was recorded with 
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. In 

(b), 



the laser set to an early time Tl of the decelerated atom 
pulse, both ^°Ne and ^^Ne are observed, but the fraction 
of ^^Ne (22 %) exceeds the natural abundance of 9.3 %. If 
the laser is set to a late time T2 of the decelerated atom 

only ^°Ne is observed 



pulse (see FIG. 15 (d)), however 
in the time-of-flight trace. 

The experimental results obtained with both decelera- 
tion pulse sequences can be understood from the phase- 
space diagrams depicted in panels (e) and (f) of FIG. 15 
which were obtained from particle-trajectory simulations 
Whereas efficient separation is possible with deceleration 
pulse sequences optimized for the lighter isotope, only 
partial separation can be achieved with the sequence op- 
timized for the heavier isotope. This overall behavior 
shows similarities to the behavior observed by Bethlem 
et al. in the multistage Stark deceleration of ammonia iso- 
topomers |25j . For a given deceleration pulse sequence, 
the particle with the smaller magnetic-moment-to-mass 
ratio (in this case ^^Ne) must lose a larger amount of 
kinetic energy to stay synchronous with the switching of 
the flelds. Therefore, it experiences a pulse sequence with 
a larger phase angle during propagation through the de- 
celerator, which leads to a smaller phase-space aceptance. 
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(c) Deceleration pulse sequence Ne 
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(b) Deceleration pulse sequence Ne 



(d) Deceleration pulse sequence Ne 
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(f) Deceleration pulse sequence Ne 
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FIG. 15. (a)-(d) Time-of-flight mass-spectrometer signal obtained by photoionizing metastable neon atoms in the beam. The 
deceleration pulse sequences used were generated for a phase angle of 40° and an initial velocity of the synchronous particle 
of 580 m/s. (e)-(f) Phase-space distributions of decelerated atoms obtained from particle-trajectory simulations, which were 
performed for an equal number of ^"Ne and ^^Ne atoms emerging from the pulsed valve. The dashed lines indicate the position 
in phase space overlapping with the ionization volume of the laser. 



Generally speaking, lighter isotopes are more strongly in- 
fluenced by the field gradients than the heavier ones and, 
unlike the heavier ones, can be decelerated even by se- 
quences optimized for different isotopes. Separation of 
the lightest isotope by multistage Zeeman deceleration is 
thus possible using deceleration sequences optimized for 
their larger magnetic-moment-to-mass ratio. 



IV. CONCLUSION AND OUTLOOK 

In this work on multistage Zeeman deceleration of 
metastable neon atoms, the properties of the decelerated 
atoms were investigated using time-of-flight and imaging 
detection techniques. These methods enabled the com- 
plete characterization of the longitudinal and transverse 



velocity distributions of the decelerated atoms. With 
their widely tunable velocity and with longitudinal and 
transverse temperatures of about 10 mK, cold atom and 
molecule samples produced by multistage Zeeman decel- 
eration are expected to find applications in the study of 
reaction dynamics at low temperatures, in precision spec- 
troscopy, and as starting point for further cooling steps. 
Multistage Zeeman deceleration leads to the production 
of samples with complete state selectivity and complete 
angular-momentum orientation as demonstrated here for 
the Mj = 2 level of metastable neon. The potential of 
such samples for reaction dynamics is evident, but so far 
unexplored. 
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